High field level crossing studies on spin dimers in the low dimensional quantum spin 
system Na 2 T2(C 2 04)3(H 2 0) 2 with T=Ni,Co,Fe,Mn 



C. Menncrich 1 , H.-H. Klauss 1 , A.U.B. Wolter 1 , S. Siillow 1 , F.J. Littcrst 1 , C. Golze 1 ' 2 , R. Klingeler 2 , 
V. Katacv 2 , B. Biichner 2 , M. Goiran 4 , H. Rakoto 4 , J.-M. Broto 4 , O. Kataeva 5 , and D.J. Price 6 
1 Institut fur Physik der Kondensierten Materie, TU Braunschweig, 
Mendelssohnstr.3, D-38106 Braunschweig, Germany 
2 Leibniz- Institute for Solid State and Materials Research IFW Dresden, P.O. Box 270116, D-01171 Dresden, Germany 

4 Laboratoire National des Champs Magnetiques Pulses, 31432 Toulouse Cedex 04, France 

5 Arbuzov Institute of Organic and Physical Chemistry, RAS, 4^0088 Kazan, Russia and 

6 WestCHEM, Department of Chemistry, University of Glasgow, Glasgow, G12 8QQ, UK 

(Dated: 6th February 2008) 

In this paper we demonstrate the application of high magnetic fields to study the magnetic 
properties of low dimensional spin systems. We present a case study on the series of 2-leg spin-ladder 
compounds Na2T2(C204)a(H20)2 with T = Ni, Co, Fe and Mn. In all compounds the transition 
metal is in the T 2+ high spin configuation. The localized spin varies from S=l to 3/2, 2 and 5/2 
within this series. The magnetic properties were examined experimentally by magnetic susceptibility, 
pulsed high field magnetization and specific heat measurements. The data are analysed using a spin 
hamiltonian description. Although the transition metal ions form structurally a 2-leg ladder, an 
isolated dimer model consistently describes the observations very well. This behaviour can be 
understood in terms of the different coordination and superexchange angles of the oxalate ligands 
along the rungs and legs of the 2-leg spin ladder. All compounds exhibit magnetic field driven ground 
state changes which at very low temperatures lead to a multistep behaviour in the magnetization 
curves. In the Co and Fe compounds a strong axial anisotropy induced by the orbital magnetism 
leads to a nearly degenerate ground state and a strongly reduced critical field. We find a monotonous 
decrease of the intradimer magnetic exchange if the spin quantum number is increased. 

PACS numbers: 75.50.Xx,75.10.Pq, 76.30.-v, 75.30.Gw 



INTRODUCTION 

In recent years, the physical properties of low dimen- 
sional spin systems have attracted a lot of attention. For 
isotropic magnetic interactions described in the Heisen- 
berg model and low spin quantum numbers the ground 
state properties are strongly influenced by quantum fluc- 
tuations and often pure quantum ground states are found 
in a macroscopic system [l|, 0, H, 13, 0] • Therefore these 
systems are ideal model systems for quantum mechan- 
ics. Low dimensional spin systems can be realized by 
linking transition metal via via organic molecules. Using 
the rich variety of organic ligands on the transition metal 
complex the dimensionality and strength of the magnetic 
interaction can be controlled. 

Within this Springer series of lecture notes in physics a 
recent volume is devoted to high magnetic field studies in 
physics In that book several contributions describe 
the physics of one dimensional magnets in high magnetic 
fields [3, Q and electron spin resonance (ESR) on molec- 
ular magnets In this article we demonstrate the ap- 
plication of high magnetic fields to study the magnetic 
properties of spin dimers and determine the parameters 
describing the system. The magnetic properties are ex- 
amined experimentally by magnetic susceptibility, pulsed 
high field magnetization, high field ESR and specific heat 
measurements. 

We present a case study on a series of structural 2-leg 



spin ladders with different transition metal ions, namely 
Na2T2(C204)3(H20)2- This series consists of four 
isostructural compounds with magnetic ions T=Ni(II), 
Co(II), Fe(II) and Mn(II) with spin 1, 3/2, 2 and 5/2 
respectively. One aim of this work is to elucidate the ef- 
fect of a gradual increase of the spin multiplicity towards 
classical magnetism. In particular, the characteristic sig- 
natures of a strong single ion anisotropy in high field 
magnetization measurements arc shown. We also discuss 
the dependence of the magnetic superexchange interac- 
tion strength on the number of 3d electrons, the influence 
of orbital moments and the validity of a description in 
the spin Hamiltonian model using isotropic Heiscnbcrg 
exchange. 



In the first and second part, we describe the chem- 
ical synthesis and crystal structure of the samples and 
give a short overview of the spin hamiltonian used in the 
analysis. In the following four parts, we discuss the dif- 
ferent compounds starting with the Ni(II) and Mn(II) 
compounds which have an orbital singlet ground state, 
followed by the Co(II) and Fe(II) compounds with an or- 
bitally degenerate ground state. Finally, we compare the 
results obtained on the four samples. 



SYNTHESIS AND CRYSTAL STRUCTURE 

The compounds NasTa^CUM^O^ with T = Ni, 
Co, Fe and Mn (we abbreviate the chemical structure 
with STOX for the general structure and SNOX, SCOX, 
SIOX and SMOX for the individual compounds) are syn- 
thesized in a hydrothermal reaction from solutions con- 
taining a very high sodium halide concentration 1^, 11 1 • 
They occur as green crystals of SNOX, purple crystals 
of SCOX, yellow crystals of SIOX and white crystals 
of SMOX. Phase homogeneity and purity of all com- 
pounds were established by a combination of optical mi- 
croscopy, powder X-ray diffraction and elemental analy- 
sis. The crystal size of the different compounds is micro- 
crystalline for the samples SCOX, SIOX and SMOX. 
Only for the SNOX compound single crystals up to 2 mg 
were grown in lower yielding reactions by slowing down 
the cooling rate, and reducing the concentration of Ni 2+ 
and (C 2 4 ) 2 -. 

Single crystal X-ray structure determination reveals 
STOX to crystallize in the monoclinic space group P2i/c 
(#14) with slightly different crystallographic axes lengths 
below 2%. The key features of the structure are the fol- 
lowing: The T(II) ion experiences a pseudo-octahedral 
coordination environment. It is coordinated in a cis ge- 
ometry by two chelating and crystallographically inde- 
pendent oxalate dianions. The remaining coordination 
sites arc filled by a monodentate oxalate oxygen atom and 
a water molecule. T-0 bond lengths lie within the range 
2.0 to 2.1 A. We note that opposite pairs of O have simi- 
lar T-0 lengths, and that the O- • • O separations between 
opposite pairs, for which the bond vectors are nearly di- 
rected along the crystallographic axes, are in the range 
of 4 A with a slightly bigger separation along the a-axes 
(fig. 1). The overall structure can be regarded as an 
anionic [T2(C204)3] 2 l ra_ network with a 1-D ladder like 
topology (fig. 1), where the metal ions form the vertices 
and the bridging oxalate ions form the bonds. The two 
oxalate ions bridge 3d ions in quite different ways. One 
oxalate forms a symmetric bis-chelating bridging mode 
(with a crystallographic inversion located at the C-C cen- 
troid) and links pairs of nickel ions (dT---T(rung) ~ 5.3 A) 
forming the rungs of the ladder. The second oxalate is 
unsymmetrically coordinated, chelating to one T(II) ion 
and forming a monodentate coordination to a second T 
ion, with 1,3-syn anti geometry through the bridging car- 
boxylate, this mode forms the ladder legs (&T-~T{leg) ~ 
5.8 A). These structural linkages are likely to provide the 
only significant pathways for magnetic superexchange. 

The low dimensionality of the dominant magnetic ex- 
change interaction in all compounds of this series is al- 
ready evident from the magnetic susceptibility measured 
on powder samples between 2 K and 300 K (fig. [2]). 
All compounds exhibit a Curie- Weiss behaviour at high 
temperatures, a maximum between as 10 K and 50 K 





Figure 1: Crystal structure of Na2T 2 (C204)3(H 2 0)2. For 
clarity, Na and H2 are not shown. 



and a strong decrease towards lower temperatures. This 
behaviour is typical for a system with a dominant anti- 
ferromagnetic interaction in less than three dimensions 
since no indications of static long range order are found 
down to 4 K. A Curie- Weiss analysis of the high temper- 
ature behaviour above 100 K yields effective magnetic 
moments close to the free ion high spin state for the Mn 
system only. While the data indicate moderate orbital 
contributions for the Ni and Fe systems, these contribu- 
tions are strong for the Co system. 




Mn ->n(Mn) =6.06 (5.92) 

Fe ->n(Fe) =5.49 (4.90) 

Co ->n(Co) =5.32 (3.87) 

Ni ->n(Ni) =3.33 (2.83) 
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Figure 2: Temperature dependence of the magnetic suscepti- 
bility of powder samples of STOX . The deduced high temper- 
ature magnetic moments (in units of fis in a Curie analysis) 
are shown in comparison with the spin only value calculated 
for the free ions (in brackets). 
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THEORETICAL MODEL 

To describe the magnetic properties we use the spin 
Hamiltonian approach. This approach neglects an ex- 
plicit description of the orbital degrees of freedom using 
pure spin coordinates. The orbital contributions treated 
as a perturbation lead to an anisotropic g-factor and 
an anisotropy of the spin orientation described by an 
anisotropy tensor D_- The spin Hamiltonian approach 
is often valid for transition metal ions because their or- 
bital moments are known to be quenched [HI, HH, Q, EH • 
It is usually working well for ions with orbitally non- 
degenerate ground states. 

Taking into account the two different exchange path- 
way topologies, we assume a stronger magnetic interac- 
tion J along the rungs than the magnetic interaction K 
along the legs. This leads to an isolated dimer approxi- 
mation where each dimer consists of two 5=1, 3/2, 2 or 
5/2 spins for Ni(II), Co(II), Fc(II) or Mn(II) respectively 
on the rungs of the ladder. Since all T ions on a ladder 
are crystallographically equivalent they share the same 
strength and orientation of the single ion anisotropy D. 
Including the Zeeman energy in the external magnetic 
field, the spin Hamiltonian of the system in the dimer 
approximation is given by 



H 



JS 1 S 2 +3/xsB(S 1 + S 2 )+ S iR S i- 



i=l,2 



Here, we assume that the magnetic field is applied 
along the z-axis, the main axis of the crystal field 
anisotropy tensor D. For pure axial anisotropy this ten- 
sor is diagonal, defined as D xx = - 1/3 D, D yy = - 1/3 D 
and D zz = + 2/3 D. The effects of an in-plane anisotropy 
could be considered in the anisotropy tensor by the pa- 
rameter E with E = \{D' XX - D' yy ). Usually \D\ > 3E 
holds. Therefore in our study we neglect the effect of 
a possible in-plane anisotropy. To calculate the energy 
levels for applied fields along different directions, we ro- 
tate the axis of the axial anisotropy tensor by an angle 7 
about a perpendicular axis. Then the spin Hamiltonian 
can be expressed as 

H = JS 1 S 2 +g yf i B B{S lz +S 2z ) + ]T S^DU^Sj (2) 

i=l,2 

with the rotation matrix U and an angle dependent 
g-factor g 7 . 

The magnetization in units of /j,b per dimer can then 
be calculated by using the equation 



dF 



oz 



M dmer = ^ = gilB o§- (3) 
e^/ T * 4 ( J, D, g, j)Sz*i(J, D, g, 7) 



where F = —ksTlnZ is the free energy, Z is the par- 
tition function and ^i(J,d, 5,7) and E{ are the eigen- 
vectors and eigenvalues of equation ©. The total mag- 
netization M tot then consists of the main contribution 
(U)) and comprises also a temperature independent term 
AIq = Mdia + M vv which includes a diamagnetic contri- 
bution Mdia and a Van-Vleck paramagnetic susceptibility 
M vv , as well as a Curie contribution C/T with a Curie 
constant C owing to paramagnetic impurities: 



Mtot = M dlmer + C/T + M . 



(5) 



To analyse the magnetization data, we developed a fit 
routine which numerically diagonalizes the Hamiltonian 
([2]) , calculates the magnetization M tot using equation ([5]) 
and varies the parameters in the Hamiltonian to minimize 
the mean square deviation between the data and this 
model. 

The isolated dimer approximation may be improved 
by introducing an effective exchange interaction K along 
the legs treated in mean field approximation in the cal- 
culation of the magnetic susceptibility x : 



Xladder — \diraer I iX ~l~ K-Xdimer)- 



(6) 



Furthermore we calculate the magnetic specific heat by 
calculating the partition function Z for the Hamiltonian 



c P . m a g = Td/dT(lnZ - 1/T (d(lnZ) / 8(3)) = (7) 
= T ^ HEeWl .I»%^)) (8) 



T 



8(3 



(4) 



with [3 = 1/fcgT, leading to an analytical function 
for a given angle 7 which can be used to analyse the 
experimental specific heat with a standard x 2 fit routine. 



Na 2 Ni2(C204)3(H 2 0)2 

Antiferromagnetic S=l chain systems are of particu- 
lar interest in quantum magnetism since they show a 
non-magnetic singlet ground state with a spin excita- 
tion gap (Haldane chains) [H, 17 1. The physical prop- 
erties of two coupled S=l chains have been studied the- 
oretically [H, 0, [2^, [2l[ but no experimental realiza- 
tion of a Haldane ladder system has been identified to 
date. Therefore we start our investigations with the 
Na2Ni2(C 2 04)3(H 2 0)2 (SNOX) compound. In this com- 
pound, the Ni(II) ions are in the 3d 8 configuration with 
an orbitally non-degenerate ground state and a total spin 
of S = 1 per ion. 
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Magnetic Susceptibility 

We examined several single crystals in magnetic sus- 
ceptibility measurements using a Quantum Design Mag- 
netic Properties Measuring System (MPMS) in external 
fields of B ex = 2 T and 5 T in the temperature range 
2-300 K. In these experiments the external field was ori- 
ented along the a-, b- and c-axes as well as along different 
intermediate angles with respect to the a-axis. 

Fig. [3] shows the temperature dependence of the mag- 
netization along the a-axis and at different angles with 
respect to the a-axis in an external field of B ex = 5 T. 
Similar to the powder measurement presented in fig. H\ 
and the single crystal measurements in B ex = 2 T 
one can clearly see a pronounced maximum for all direc- 
tions with a strong downturn below « 50 K suggesting a 
non-magnetic spin-singlet ground state. This behavior in 
general is expected as well for an isotropic two-leg spin 
ladder as for the two limiting cases, an isolated S = 1 
Haldane chain or a system of antiferromagnctically cou- 
pled dimers 0, M, M, EI H E | . 

The inset of fig. [3] shows the behavior at very low 
temperatures where the magnetization along a is smaller 
and approaches zero for T — > in contrast to directions 
transverse to a. Since the a-axis is close to the axial 
anisotropy axis as we will see below this indicates an 
easy axis of the Ni spin moments. The difference between 
experiment and theory at very low temperatures is due 
to an overestimated temperature independent Van Vleck 
contribution in the theoretical fit. 
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Figure 3: Temperature dependence of the magnetization of a 
SNOX single crystal at an external field of B ea! = 5 T along 
the a-axis (circles) and at different angles with respect to the 
a-axis (triangles and squares). The solid line represents the 
fit as described in the text. 

To analyse the magnetization data, we performed a 
combined fit of Mtot for B ex = 5 T along all mea- 
sured directions using an isotropic magnetic exchange 
coupling constant J and a fixed absolute value of the 
axial anisotropy of \D\ = 11.5 K. This anisotropy value 



was determined by the S=l zero field splitting measured 
by ESR as described below. From these measurements 
we obtain an intradimer coupling constant of J = 44 K, 
an interdimer coupling constant of K = K and ^-values 
of g a = 2.215, 545 = 2.300 and g 90 = 2.330 . The fitted 
angles to the anisotropy axis are j a = 19°, 745 = 57° 
and 790 = 90°. In addition a temperature independent 
contribution \o = 0.03 /ig/dimer has been used. The 
fit results are shown in fig. [3] as solid lines. They are 
in good agreement with the results obtained in llj for 
B ex — 2 T. 

To determine the absolute value of the axial anisotropy 
constant D we have performed tunable high-field ESR 
measurements of SNOX at frequencies up to v ~ 1 THz 
in magnetic fields up to B ~ 40 T. Details of experimen- 
tal set-ups can be found in Ref . [22[ ■ A complex spectrum 
comprising a main line and a number of weak satellites 
has been observed. The ESR-intensity shows a thermally 
activated behavior similar to the T-dependence of the 
static magnetization (see fig. 01 inset) thus ensuring that 
the ESR response is determined by the bulk Ni-spins. 
The frequency vs. magnetic field dependence of the ESR 
modes is shown in the main panel of fig.[H In our analysis 
of the ESR data we focus the attention on the strongest 
resonance line 1 which i/(B)-dependence exhibits the in- 
tercept with the frequency axis at vq » 239 GHz (details 
of the ESR analysis can be found in Ref. [HI). This inter- 
cept implies a finite energy gap A = voh/kg = 11. 5K 
for this resonance excitation. This zero field gap can be 
straightforwardly identified with the zero field splitting of 
the first excited 5=1 triplet state of the dimer and gives 
directly the magnitude of the axial anisotropy parameter 
A = LDI = 11.5 K. 



The knowledge of j a , 7^ and 7 C [llj allows a deter- 
mination of the orientation of the anisotropy axis in the 
crystal. It is oriented parallel to the connecting line of the 
opposite nearest neighbor oxygen ions along the a-axis, 
tilted by 18° with respect to the crystallographic a-axis 
mainly towards the b axis [ll||. Note, that analyzing the 
data with a negative D =-11.5 K leads to a similar good 
agreement between the model and the data. However, 
the positive D value is strongly supported by the high 
field magnetization data described in the next section. 

The relative energies of the spin states of SNOX calcu- 
lated in the framework of the Hamiltonian @ with the 
parameters yielding the best fit to the temperature de- 
pendent magnetization are plotted in fig. [5] for magnetic 
fields oriented parallel and perpendicular to the main axis 
of the anisotropy tensor. One can see that at zero field 
the spin singlet ground state S = is well separated 
from the S = 1 triplet and the S = 2 quintet state. 
The separation energy between singlet and triplet corre- 
sponds to the coupling constant J, whereas that between 
singlet and quintet corresponds to 3J. The splitting of 
the excited S — 1 and S = 2 states in zero magnetic 
field is induced by the anisotropy D. Both triplet and 
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Figure 4: Main panel: Frequency v vs. resonance field 
B res dependences (branches) of the strongest line 1 and weak 
satellites 2,3 and 4 comprising the ESR spectrum of SNOX. 
Note the intercept of branch 1 with the frequency axis which 
yields the magnitude of the anisotropy parameter D. Solid 
fine is a theoretical fit according to Hamiltonian (2) with 
\D\ — 11.5 K. Dashed lines are a linear approximation of weak 
satellite branches. Grey bar indicates the spread of the over- 
lapping satellite resonances in the high field regime. Inset: 
T-dependence of the strongest ESR line 1. Note its similarity 
to the T-dependence of the static magnetization. 



quintet levels split with increasing magnetic field due to 
the Zeeman effect. The interplay of the zero field split- 
ting and the Zeeman splitting in the determination of the 
high field properties of SNOX will be discussed in detail 
in the next section. 



High field magnetization 

As can be seen from the Brcit-Rabi diagram of the 
spin states of SNOX (fig. the splitting of the states 
in a magnetic field yields a level crossing of the ground 
state with the lowest triplet state | — ) at a field Bci 
and a second level crossing of the | — } state with the 

| } quintet state at field Bci with Bci < Bci- This 

leads to a step-like behavior in high field magnetization 
measurements. If the field is applied along the z axis, 
the critical field Bci can be estimated using the equa- 
tion Bci = (J — D/3)/(g/Xs/fcs), depending strongly 
on the sign of D. For a field perpendicular to the local 
anisotropy axis, a negative D (leading to higher criti- 
cal fields for B || z) pushes the critical field to lower 
fields and vice versa for a positive D. Since in powder 
measurements the perpendicular situation dominates the 
spectrum of the spin states, the field dependent magne- 
tization, among others, can determine the sign of the 
anisotropy D. 

The critical field strengths Bci depend on the mag- 
netic exchange J. An exchange energy of J=44 K cor- 
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Figure 5: Relative energies of the spin states, calculated for 
the magnetic field applied parallel and perpendicular to the 
z axis of the uniaxial anisotropy tensor: the triplet S — 1 
and the quintet S = 2 are well separated from the S = 
ground state by an activation energy of | J| = 44 K and 3 | J|, 
respectively. Both multiplets exhibit a zero field splitting due 
to crystal field anisotropy. The shaded areas mark the field 
ranges of ground state level crossings around 29 and 60 T in 
powder magnetization measurements (for details see text). 



responds roughly to Bci =30 T for Bci- Therefore high 
magnetic fields are needed to observe the magnetization 
steps. We performed measurements at several tempera- 
tures (1.47 K, 4.2 K, 10 K and 26 K) on a powder sample 
with a mass of 23.2 mg in magnetic fields up to 55 Tesla. 
The results are shown in fig. [5J In the low tempera- 
ture experiments the first magnetization step is clearly 
seen, whereas it smoothes at higher temperatures due 
to the thermal averaging process. The critical field Bci 
can be obtained from the derivation of the magnetiza- 
tion dM/dB for T = 1.47 K (see inset of the upper panel 
of fig. EH, T = 4.2 K and T = 10 K. At all tempera- 
tures a sharp peak is observed at 29 Tesla. Additionally 
a second increase of the magnetization at much higher 
fields of about 50 Tesla can be anticipated. This step 
is clearly observed at T — 1.47 K and T = 4.2 K and 
strongly broadened in the 10 K measurement. At 26 K 
the magnetization steps are not observed and a linear 
field dependence of the magnetization is found. 

The solid lines in figE] describe simulations using the 
dimer model with the calculated g values <?j| = 2.201 
(see above) and g± = g^Q. For the simulation of a pow- 
der measurement in this case a weighted average over 
a full set of angles between the anisotopy axis and the 
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field direction from 0° to 90° in steps of 1° is performed 
with varying g from g\\ to g±. Since high magnetic fields 
lead to strong nonlinearities in the magnetization curve 
this approach is more realistic than the two step average 
M = (My+2 M]_)/3. For a good description of the high 
field magnetization we need to use a coupling constant 
of J = 44.5 K instead of J = 44 K. A small deviation is 
found in the field range of 30 - 33 T in the T = 1.47 K 
measurement, which may be caused by a little shift of 
the sample induced by the external field. 

The simulations for a negative D = -11.5 K did not 
describe the data reasonably well for any J. In this case 
the simulated magnetization step is clearly broader than 
the step of the experimental data set. The enhanced 
broadening of the magnetization step for D = -11.5 K 
can be explained by the resulting switching of the m = 
state with the |m|=l states which leads to different level 
crossing field strengths Bci- 
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Figure 6: High field magnetization of a SNOX powder sample. 
Upper panel: Measurements at 1.47 K and 4.2 K, the inset 
shows the differential magnetization dM/dB at 1.47 K. Lower 
panel: measurements at 10 K and 26 K, the inset shows the 
simulations for T = 0.01K up to 80 T field strength. The 
solid line represents the simulation as described in the text. 



B C i = 29 T. In our simulations, we get B C i = 27.8 T 
for B parallel to the magnetic anisotropy axis and 
Bci = 30.2 T for B perpendicular to the magnetic 
anisotropy axis, resulting also in Bci = 29 T in a pow- 
der average. Notice, that the powder averaging leads to 
a broadening of the magnetization steps. For the sec- 
ond step, we get B C i = 67.3 T and B C i = 51.1 T for 
B parallel and perpendicular, respectively. This differ- 
ence leads in a powder sample at the lowest temperature 
of our experiment, T = 1.47 K, to a continuous dou- 
ble sigmoid-like increase of the magnetization in the field 
range of 51 - 67 Tcsla, of which only the first increase is 
observed in the measurement. 



Specific Heat 

As shown in fig. [5] the energy gap A between the sin- 
glet ground-state and the | — ) state of the magnetic S = 1 
triplet is reduced by an external field from Aqt ~ 41.1 K 
to zero for a field perpendicular to the anisotropy axis 
at a level crossing field of Bci ~ 30.2 T. Then at higher 
field the gap increases with further increasing field un- 
til the | ) state of the S = 2 triplet anti-crosses the 

S = singlet at « 42 T . At higher fields the gap de- 
creases again towards zero at the second ground state 
change at Bc2 ~ 51.3 T. This behaviour can be directly 
observed in the calculated magnetic specific heat accord- 
ing to equation ©(fig. \7§. At low fields only one broad 
maximum is found. Increasing the field leads to a de- 
crease of the gap. Approaching Bci a low temperature 
peak is separated from the main component. At the crit- 
ical field this low temperature peak sharpens and shifts 
towards zero temperature whereas the maximum of the 
main component shifts to slightly higher temperatures. 
For fields above Bci the low temperature peak shifts to 
higher temperatures. 




temperature (K) 



From the differential magnetization dM/dB we can 
determine the critical field Bci of the powder sample to 



Figure 7: Calculated temperature dependence of the magnetic 
specific heat of SNOX for high external fields. 
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Measurements of the specific heat c p (T) were per- 
formed on a single crystal of mass 2.04 mg using a Quan- 
tum Design PPMS in zero-field and an external field of 
B = 9 T perpendicular to the a-axis in the temperature 
range 2 - 300 K (fig®. The measurements differ only in 
the low temperature regime where the magnetic specific 
heat is dominant 1231 . 



oxalate-bridged dimers 0, [H, [H, [27|, [H, 0, H 
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Figure 8: Temperature dependence of the specific heat of a 
SNOX single crystal (2.04 mg) at zero field (open circles) 
and B ea; = 9 T (filled squares) compared with the calculated 
magnetic specific heat for zero field (thick line) and B ex = 9 T 
(thin line). The inset shows the measurements on an extended 
temperature range. 

Due to the large coupling constant of the dimer, the 
magnetic contribution to c p (T) is strongest in the tem- 
perature range around T « 15 - 50 K. Since the De- 
bye temperature Tp of such metal-organic compounds 
is typically in the range of fOO - 200 K, the magnetic 
contribution to c p cannot be extracted with a standard 
T 3 -fit for the lattice contribution which is possible only 
for T <C T/j. By analyzing the difference of both mea- 
surements c p (9T) - c p (0T) as shown in figJ5]we subtract 
the phonon contribution. A fit using the isolated dimer 
model with variable parameters J and g in the temper- 
ature range T < 19 K results in a very good agreement 
for a coupling constant J = 43.7(5) K and a g-value of 
2.41(5) (2|. 

In conclusion, the specific heat measurements on 
SNOX are consistent with the susceptibility and high 
field magnetization measurements. All experiments show 
that the isolated dimer model in the spin hamiltonian 
approach provides a very good description of the sys- 
tem. This can be understood considering the different 
coordination and superexchange angles of the exchange 
mediating oxalate molecule along the rungs and legs of 
the ladder. A detailed discussion is presented below. The 
derived parameter values for D, J and the g- factor are 
typical for Ni(II) in a pseudo-octahedral environment and 
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Figure 9: Specific heat difference c p (9T) - c p (0T) of SNOX. 
The inset shows the difference in an extended temperature 
range. 



Na 2 Mn2(C204)3(H 2 0)2 
Magnetic Susceptibility 

InNa 2 Mn 2 (C 2 04)3(H20)2 the Mn(II) (3d 5 ) configura- 
tion results in an orbitally non-degenerate ground state 
and a total spin of S = 5/2 per ion. This leads to a very 
specific behaviour compared with the Ni(II), Co(II) and 
Fc(II) ions: Whereas for the ions with a more than half- 
filled 3d shell the orbital magnetic moments are quenched 
and have to be considered as a perturbation, the orbital 
angular momentum L for the 3d 5 configuration is zero 
following Hund's rules. Therefore we can neglect the 
anisotropy term and the Hamiltonian simplifies to 



H = JS 1 S 2 + 5 /x B B(Si + S 2 ). 



(9) 



We measured a powder sample of mass 24.8 mg in a 
home built vibrating sample magnetometer 3l| in exter- 
nal fields of B ex = 1 T and 16.8 T in the temperature 
range 2-300 K. The measurements are shown in fig. lfTO")) . 
For B ex = 1 T a clear maximum is observed at around 
10 K whereas for B ex = 16.8 T the magnetization in- 
creases monotonously with decreasing temperature indi- 
cating that at this field strength the antiferromagnetic 
intradimer interactions are suppressed by the strong Zee- 
man interaction of the individual spins in the external 
field. 

To analyse the magnetization data, we performed a 
combined fit of M tot for B ex = 1 T and B ex = 16.8 T 
using an isotropic magnetic exchange coupling constant 
J and an averaged g-factor. Since without anisotropy 
the magnetizations for applied external fields parallel and 
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Figure 10: Temperature dependence of the magnetization of a 
SMOX powder sample at external fields of B ex i = 1 T (trian- 
gles) and B ea; i = 16.8 T (squares). The 16.8 T measurement 
is divided by the field. The solid line represents the fit as 
described in the text. 



perpendicular to the z-axis are identical we performed an 
analysis without powder averaging. We obtained an in- 
tradimer coupling constant of J = 3.5 K and a g- value 
of g = 2.01. The fit is done without a temperature in- 
dependent contribution and without considering param- 
agnetic impurities. The results are shown in fig. 1101 as 
solid lines. To prove the validity of the isotropic model 
(i.e. the & S§/2 ground state) we performed an analysis of 
the B ex = 1 T measurement including an anisotropy D 
and an anisotopic g-factor in a powder average. This ap- 
proach did not improve the quality of the fit. The same 
holds for including an interdimer coupling K on the mean 
field level. 

The different behaviour of the 1 T and the 16.8 T 
magnetization curves can be explained considering the 
energy levels of the system. The Breit-Rabi diagram of 
SMOX shows 6 multiplets, one for each spin between 
S = and S = 5. These multiplets are degenerate for 
B = reflecting the absent anisotropy D. In an exter- 
nal magnetic field the multiplets split due to the Zeeman 
interaction and five ground state level crossings are ob- 
served: for J = 3.5 K and g = 2.01 at approximately 2.6, 
5.2, 7.8, 10.4 and 13 T. A non-magnetic ground state 
(i.e. dE/dB = 0) is only found for fields below 2.6 T. 
At higher fields a magnetic ground state exists with a 
finite magnetization for T — > which increases steplikc 
at the level crossings. The different ground states can be 
observed in the low temperature behaviour of the temper- 
ature dependent magnetization curves. The temperature 
dependence of the magnetization for fields below 2.5 T 
shows a strong downturn to zero reflecting the depopu- 
lation of the magnetic levels whereas for higher fields a 
finite value of 2, 4, 6, 8 and 10 [Ib per dimcr is found 
for the different field ranges. For the highest field regime 



above 13 T, the magnetization increases monotonously 
with decreasing temperature to the maximum value of 
10 fiB per dimer. 



High field magnetization 

Field dependent magnetization measurements were 
also performed in the vibrating sample magnetometer on 
a powder sample at temperatures 2.4, 6, 10 and 30 K 
in fields up to 16.8 T. The measurements are shown in 

fig. CD. 

The analysis is done using the Hamiltonian ((9]). The 
eigenvalues can be calculated directly using 



E,, 



J/2 (S(S+i)_(5 1 (5i+l)-S , 2 (S , 2+l))- 



■g mu B B {S lz +S 2z ) 
(10) 



were the basis | S m s ) is given by total spin 
S = Si + S2 = 0... 5 with magnetic quantum numbers 
m s from -S to S leading to 36 eigenvalues. 

This basis allows to calculate the exact eigenvalues 
with parameters J, g and B leading to an analytical func- 
tion for the magnetization -Mdj mer (J, g, B) using equation 
©. With this function Md imer (J, g, B) we performed a 
combined fit for all temperatures leading to a coupling 
constant J = 3.5 K and g = 2.04 in excellent agreement 
with the temperature dependent measurements. The 
fit results are shown in figfTT] as solid lines. The grey 
solid line shows the calculated magnetization for a very 
low temperature T = 0.1 K. Only at this temperature 
the magnetization steps at the level crossing fields are 
present. These steps are smoothed even at the lowest 
measured temperature of 2.4 K due to the thermal aver- 
aging process and the small coupling constant J. 

The deduced parameter value for J and the g-factor are 
typical for Mn(II) in a pseudo-octahedral environment 
and oxalate- bridged dimers 32, 33, 34, 35, 36j |. 



Na 2 C02(C2O 4 )3(H 2 O)2 
Magnetic Susceptibility 

The Co( II) compound has been synthesized by Price 
ct al. [l(J. In SCOX, the Co(II) ions are in the 3d 7 
configuration with an orbitally degenerate ground state 
Tig and a spin of S = 3/2 per ion in the high-spin state. 
The presence of an additional orbital moment is already 
evident from the high temperature Curie- Weiss analysis 
presented in the introduction. Therefore, a description of 
the magnetic properties using the spin hamiltonian ap- 
proach will result in an anisotropic g-factor very different 
from g=2 and a strong uniaxial anisotropy due to spin- 
orbit coupling. Alternatively interacting Co(II) ions in 
octahedral environments are described by a pseudo spin 
1/2 and a strong Ising type anisotropy in the magnetic 
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Figure 11: High field magnetization of a SMOX powder sam- 
ple. The figure shows measurements at 2.4 K, 6 K, 10 K and 
30 K. The black solid lines represents the fit as described in 
the text and the grey solid line is a calculated magnetization 
for T = 0.1 K. 
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Figure 12: Temperature dependence of the magnetization of 
a SCOX single crystal at an external field of T5 ex — 0.01 T 
along the a-axis (circles) and perpendicular to the a-axis (di- 
amonds). The data are taken from ref. [37j], the solid lines 
represent an isolated dimer fit as described in the text. 



interaction. Since the purpose of this work is a compar- 
ison of the different spin multiplicities on the transition 
metal site wc use the model hamiltonian equation 

The first magnetic susceptibility data on a single crys- 
tal parallel and perpendicular to the a-axis have been 
performed by Honda et al. [37j |. There the data have 
been analyzed using a phcnomenological S=l/2 two-leg 
spin ladder function in a limited temperature range only 
(7 to 20 K). From our systematic study of STOX it is 
evident that also the Co(II) system can be described as 
nearly uncoupled dimcrs of S=3/2 spins on the rungs of 
the ladder. Therefore we present a new analysis of the 
single crystal magnetic susceptibility data based on this 
model similar to the Ni(II) system using hamiltonian (J5J). 

A combined fit along both directions (solid lines in 
fig. fT2|) results in an intradimcr coupling constant of 
J = 13.8 K, an anisotropy of D = -6.3 K and ^-values 
of g a = 2.60, and g± = 2.27 . For comparison with 
these single crystal measurements, we performed a sus- 
ceptibility measurement at B = 1 T on the powder sam- 
ple used for the high field magnetization measurements. 
An analysis using equation ([5]) with a powder average 
(M = (M||+2 M_l)/3) results in J = 11.8 K and D = -9 K 
and a large g- value anisotropy with ^-values of = 3.90, 
and g± = 1.50 . Note that J is smaller by a factor of « 
4 than in the Ni(II) system. The variation of the g-value 
is opposite to SNOX due to a different orbital contribu- 
tion to the total magnetic moment. Since no high field 
high frequency ESR measurements have been performed 
so far, the single ion anisotropy strength D had to be 
used as a free parameter. A satisfactory description of 
the susceptibility and high field magnetization data is 
obtained only for a negative value of D. In view of the 
limited accuracy of the single crystal fit the angles to the 



anisotropy axis j a = 28° and 7j_ = 38° arc consistent 
with an orientation of the anisotropy axis as in SNOX 
parallel to the connecting line of the opposite nearest 
neighbor oxygen ions along the a-axis. 

Fig. [T3l shows the low energy spin states of SCOX cal- 
culated for magnetic fields oriented perpendicular to the 
main axis of the anisotropy tensor. This orientation gives 
the dominant contribution in a powder magnetization 
measurement. For this plot we used the set of param- 
eters J— 11.8 K, D—-9 K and g±=l. 5 derived from the 
fit of the powder measurements which is also used for 
the simulation of the high field magnetization measure- 
ment as described below. Note the near degeneracy of 
the |m=0> state of the 5=1 triplet and the 5 = 
singlet due to the particular combination of J and D. 
Therefore a series of ground state level crossings is pre- 
dicted at critical field strengths Bet of « 11.6 T, 26.6 T 
and 41.7 T in this field orientation. 



High field magnetization 

We performed a pulsed high magnetic field magnetiza- 
tion measurement at 4.2 K on a powder sample of SCOX 
in magnetic fields up to 52 Tesla (fig [Tij) . It shows one 
broadened magnetization step from w 0.5 [Ib to nearly 
5.5 Hb between 10 T and 20 T in addition to a strong 
finite magnetic susceptibility (gradient in M(B)) at low 
fields. This behavior is in stark contrast to a three step 
function which can be expected for a dimer of S=3/2 
spins from the Breit-Rabi diagram (fig. [T3|) . In this 
case when the antiferromagnetic exchange is of similar 
strength as the single-ion anisotropy it is necessary for 
a qualitative and quantitative description of a powder 
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Figure 13: Low energy spin states in SCOX for a mag- 
netic field applied perpendicular to the z axis of the uniaxial 
anisotropy tensor: the |m=0> state of the S = 1 triplet and 
the 5 = singlet are nearly degenerate due to the particu- 
lar combination of J and D. The arrows indicate the ground 
state level crossing fields. 
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Figure 14: Pulsed high field magnetization of a powder sample 
of SCOX at 4.2 K. The solid lines represent the simulations 
at 4.2 and 0.1 K. The inset shows single crystal simulations 
for three different angles between the magnetic field and the 
anisotropy axis in comparison with the powder average for 
T=0.1 K. 



magnetization measurement to perform a weighted aver- 
age over a full set of angles between the anisotropy axis 
and the field direction from 0° to 90° in steps of 1° with 
varying g from g\^ to g± since the Breit-Rabi diagram de- 
pends strongly on this angle. This is illustrated in the 
inset of fig. [TJ] where we plot the calculated single crys- 
tal magnetization curves at T = 0.1 K for three different 
angles in comparison with the powder average. Due to 
the strong anisotropy for B || z the first level crossing 
at 10 T leads from S = to the |m=-3> of the S = 3 
fully polarized sextet state resulting in a sharp step of 
the magnetization curve. This situation contributes sig- 
nificantly to the powder average. The black solid line in 
fig. Q3] is calculated for T~ 4.2 K using the set of pa- 
rameters obtained by the powder susceptibility with an 
optimized value <7|| =3.0. Due to the dominance of the 
large step of the B\\z direction, the influence of the steps 
in the Biz direction in the powder sample is notice- 
able only at very low temperatures. Note that in figfTH 
at T = 0.1 K the magnetization shows two small kinks 
at the level crossing fields of the B _L z direction (arrows 
on grey solid line). These kinks result from the powder 
averaging which leads to a smoothing of the steps (cf. in- 
set of fig IT4"|) . Recently we became aware of experiments 
on Co-Oxalate by Y. Nakagiwa et al. which qualitatively 
agree with our findings [381 ] - 



In the literature a very wide range of values for D and 
the g-factor are found for Co(II) complexes. See e.g. the 



review by R. Boca [151 ] or examples in [3£ 

II SB | . 
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Na 2 Fe2(C204)3(H 2 0) 2 
Magnetic Susceptibility 

The magnetic susceptibility of a powder sample of the 
Fe( II) compound has been measured by Kreitlow et al. 
Hi]. In SIOX, the Fe(II) ions are in the 3d 6 configura- 
tion with an orbitally degenerate ground state 5 T 2g and a 
spin of S = 2 per ion in the high-spin state. The presence 
of an additional orbital moment is already evident from 
the high temperature Curie- Weiss analysis of the mag- 
netic susceptibility. Comprehensive single crystal mag- 
netic susceptibility data parallel and perpendicular to the 
a- axis have been published by Kikkawa et al. [13] show- 
ing a very strong anisotropy. Only the data parallel to 
the a-axis have been analyzed using a singlet-triplet sys- 
tem including a uniaxial anisotropy resulting in a singlet- 
triplet gap of 51 K and a very large anisotropy D=41 K. 
Here we present a new analysis of the single crystal mag- 
netic susceptibility data parallel and perpendicular to the 
a-axis based on a dimer model using hamiltonian ([2"]). 

The combined fit along both directions (solid lines in 
fig. I15p results in an intradimcr coupling constant of 
J = 6.5 K, g- values of g a = 2.67, and g± = 1.76 . Note 
that J is reduced by a factor of m 2 with respect to 
the Co(II) system. The single ion anisotropy strength 
D has been determined to D=-7.1 K. Similar to the 
SCOX compound, we performed magnetic susceptibil- 
ity measurements on a SIOUX powder sample used for 
the high field magnetization measurements and analysed 
these data in the same way. This results in a param- 
eter set of J = 6.7 K, D = - 8.7 K, g\\ = 2.67, and 
g± = 2.65 which we used to calculate the Breit-Rabi di- 
agram (see inset of fig. [TB]) . Therefore in this system the 
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Figure 15: Temperature dependence of the magnetization of 
a SIOX single crystal at an external field of B ex — 0.01 T 
along the a-axis (circles) and perpendicular to the a-axis (di- 
amonds) . The data are from ref . [43] , the solid line represents 
an isolated dimer fit as described in the text. 
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Figure 16: Pulsed high field magnetization of a powder sample 
of SIOX at 4.2 K. The solid lines represent simulations at 4.2 
and 0.1 K. The inset shows the low energy states in a Breit- 
Rabi diagram for a magnetic field applied perpendicular to 
the 2 axis of the uniaxial anisotropy tensor. 



MAGNETIC EXCHANGE PATHWAYS 



isotropic magnetic exchange and the single ion anisotropy 
are of similar size leading to the strong anisotropy in the 
measured magnetic susceptibility. In view of the lim- 
ited accuracy of the fit the angles to the anisotropy axis 
7a = 35° and "f± = 52° are again consistent with an ori- 
entation of the anisotropy axis parallel to the connecting 
line of the opposite nearest neighbor oxygen ions along 
the a-axis. 



High field magnetization 



We performed a pulsed high magnetic field magnetiza- 
tion measurement at 4.2 K on a powder sample of SIOX 
in magnetic fields up to 52 Tcsla (fig fTo)) . It shows the 
same oveall features as the experiment on the Co ana- 
logue, namely a broadened magnetization step from ss 
1.0 hb to nearly 10 p,B per dimer between 7 and 20 T in 
addition to a strong finite magnetic susceptibility at low 
fields. The low energy part of the Brcit-Rabi diagram 
for B perpendicular to the anisotropy axis using a set of 
parameters J— 6.7 K, D=-8.7 K, and g±=2.65 optimized 
to describe the high field magnetization data is shown in 
the inset of fig[l6] The |m=0> state of the £ = 1 triplet 
and the S = singlet are very close and several level 
crossings occur between 5 and 40 T. 

In the literature similar to Co (II) a very wide range 
of values for D and the (/-factor is found for Fe(II) com- 
plexes. See e.g. the review by R. Boca [l5| or examples 
in 0, H 0, 1H, 0,0,0,|55|. 



An unexpected result of the present study is the very 
large difference between the strength of the magnetic ex- 
change interaction along the rungs and legs of the spin 
ladder leading to the applicability of the dimer model in 
this system. This can be qualitatively understood consid- 
ering the different coordination and superexchange angles 
of the exchange mediating oxalate molecules. 

On the rungs the oxalate forms a ^-1,2,3,4 bridge be- 
tween two transition metal ions. This bridge provides 
two symmetric superexchange pathways. On each path 
the transition metal 3d x 2 -y 2 orbitals have an enhanced 
electron probability density extending directly towards 
the corresponding oxygen ions. There is a direct overlap 
of the transition metal 3d x 2 -y 2 and the O 2p wave func- 
tions forming a bonds. The polarized 2p orbitals them- 
selves are strongly overlapping. Therefore the interme- 
diate carbon atom is not involved in the superexchange 
mechanism resulting in a strong antiferromagnctic su- 
perexchange interaction (see fig. [T7|) . In the literature 
several /x-1,2,3,4 oxalato-bridged Ni(II) dimer and chain 
systems are reported (see e.g. [3(1 0] and references 
therein) with J values in the range of 20 to 42 K. 

For the oxalate bridge along the legs of the spin ladder 
the situation is completely different. In this case a /x-1,2,3 
oxalato bridge is formed. As shown in the right panel of 
fig. [T7]only one transition metal ion forms two covalent 
bonds (one with the x 2 -y 2 orbital and one with the 3z 2 
- r 2 orbital) with the oxalate molecule whereas of the 
second transition metal ion only the 3z 2 - r 2 orbital forms 
one covalent bond. The O 2p orbitals involved in the T-0 
bonds are not overlapping with each other and therefore 
either the carbon atoms or orthogonal O 2p orbitals are 
involved which strongly suppresses the strength of the 
superexchange mechanism and may even lead to a weak 
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Figure 17: T-(ox)-T coordination and proposed magnetic 
exchange pathways on the rungs (left panel) and on the legs 
of the spin ladder (right panel). The involved T 3d and oxygen 
2p orbitals are shown in black. 



ferromagnetic exchange. Similar /i-1,2,3 oxalato bridges 
are reported for Cu(II) systems with magnetic exchange 
strengths in the range of -0.2 to 0.3 K 57, [EH consistent 
with the estimate of \K\ [ll[ being negligibly small in 
STOX. 

A second interesting result is the continuous decrease 
of the isotropic magnetic exchange strength J going from 
3d 8 (Ni(II)) to 3d 5 (Mn(II)) shown in fig. OH The same 
systematics is present e.g. in the series of three dimen- 
sional antifcrromagncts KTF3 with T=Ni,Co, Fe, Mn 
[59l | . There the nearest neighbor exchange varies from 
44-50 K in KNiF 3 over 19 K in KC0F3 and 6 K in KFeF 3 
to 3.6 K in KMnF3. This systematics can be regarded as 
a result of the description of the total magnetic exchange 
energy in terms of the product JS1S2. As described 
above the 3d x 2 -y 2 orbital contributes the most to the 
total spin exchange energy. In Co(II), Fe(II) and Mn(II) 
the additional spin density is located in the t2 9 orbitals. 
These orbitals do not point directly towards the oxygen 
ions. They may contribute only very little to the mag- 
netic exchange energy but increase the total spin value. 
As a consequence the numerical value of the exchange 
constant J is reduced. 



CONCLUSION 

In this paper we demonstrate the application of 
high magnetic fields to study the magnetic proper- 
ties of low dimensional spin systems. We present a 
case study on the series of 2-lcg spin-ladder compounds 
Na 2 T2(C204)3(H 2 0)2 with T = Ni, Co, Fe and Mn. In 
all compounds the transition metal is in the 2+ high spin 
configuation. The localized spin varies from S=l to 3/2, 
2 and 5/2 within this series. The magnetic properties 
were examined experimentally by magnetic susceptibil- 
ity, high field ESR, pulsed high field magnetization and 
specific heat measurements. The data are analysed us- 
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Figure 18: Change of the isotropic magnetic exchange con- 
stant J with decreasing number of 3d electrons, i.e. increasing 
spin quantum number, in the series Na2T2(C204)3(H20)2- 



ing a spin hamiltonian description with Heisenberg ex- 
change interaction. Although the transition metal ions 
form structurally a 2-leg ladder, an isolated dimer model 
consistently describes the observations very well. All 
compounds exhibit magnetic field driven ground state 
changes which at very low temperatures lead to a mul- 
tistep behaviour in the magnetization curves. In the 
Co and Fe compounds a strong axial anisotropy induced 
by the orbital magnetism leads to a nearly degenerate 
ground state and a strongly reduced critical field. We 
find a monotonous decrease of the intradimcr magnetic 
exchange if the number of 3d electrons is decreased which 
indicates that the additional spin density in the t2 9 or- 
bitals does not contribute to the exchange energy. 
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